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ABSTRACT 

Exploiting the deep high-resolution imaging of all five CANDELS fields, and accurate redshift information provided 
by 3D-HST, we investigate the relation between structure and stellar populations for a mass-selected sample of 
6764 galaxies above 10 u) Mq, spanning the redshift range 0.5 < z < 2.5. For the first time, we fit two-dimensional 
models comprising a single Sersic fit and two-component (i.e., bulge + disk) decompositions not only to the //-band 
light distributions, but also to the stellar mass maps reconstructed from resolved stellar population modeling. We 
confirm that the increased bulge prominence among quiescent galaxies, as reported previously based on rest-optical 
observations, remains in place when considering the distributions of stellar mass. Moreover, we observe an increase 
of the typical Sersic index and bulge-to-total ratio (with median B/T reaching 40%-50%) among star-forming 
galaxies above 10 11 Mq. Given that quenching for these most massive systems is likely to be imminent, our 
findings suggest that significant bulge growth precedes a departure from the star-forming main sequence. We 
demonstrate that the bulge mass (and ideally knowledge of the bulge and total mass) is a more reliable predictor 
of the star-forming versus quiescent state of a galaxy than the total stellar mass. The same trends are predicted by 
the state-of-the-art, semi-analytic model by Somerville et al. In this model, bulges and black holes grow hand in 
hand through merging and/or disk instabilities, and feedback from active galactic nuclei shuts off star formation. 
Further observations will be required to pin down star formation quenching mechanisms, but our results imply that 
they must be internal to the galaxies and closely associated with bulge growth. 
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1. INTRODUCTION 

The mechanism driving the shutdown of star-forming galaxies 

(SFGs), “quenching,” remains one of the least understood 
puzzles in galaxy formation to date. In the low-redshift universe, 

galaxies show a bimodal color distribution, accompanied by 

a bimodality of morphologies (see Kauffmann et al. 2003; 

Strateva et al. 2001). Spiral galaxies have low bulge-to-disk 
ratios and are commonly the site of active star formation 

leading to blue colors. Passive galaxies are observed to be 

mostly spheroid-dominated. The color bimodality has also 
been observed at higher redshifts (e.g., Brammer et al. 2009; 
Whitaker et al. 201 1). These observations have been interpreted 

as evolutionary paths, in which one or several quenching 
processes cause the SFG to become red and passive on a 

short timescale (see Bell et al. 2004; Faber et al. 2007). 

Studies of the shape of the mass function of passive galaxies 

over cosmic time, using large surveys such as Sloan Digital 

Sky Survey (SDSS), NMBS, zCOSMOS, UltraVISTA, and 
zFOURGE, showed that the probability for a galaxy being 


quenched increases with its mass (Peng et al. 2010; Brammer 
et al. 2011; Ilbert et al. 2013; Muzzin et al. 2013; Woo et al. 
2013; Tomczak et al. 2014). Further evidence for quenching also 
comes from studies employing abundance matching techniques, 
where the cumulative abundance of galaxies is matched to that of 
halos using the results of cosmological dark matter simulations 
alongside observational constraints on the stellar mass function. 
Those infer low baryon fractions within dark matter halos, 
hinting at an efficient quenching mechanism associated with 
significant gas mass loss for both low and high mass galaxies 
(Moster et al. 2010; Behroozi et al. 2010). 

Several quenching mechanisms have been proposed, which 
act to either remove the gas from the galaxy or prevent the 
existing/inflowing gas within the galaxy to form stars. For the 
high-mass regime of galaxies, powerful active galactic nucleus 
(AGN) feedback, which may be induced by galaxy mergers (e.g., 
Hopkins et al. 2006) or internal evolutionary processes triggered 
by (violent) gravitational disk instabilities (e.g., Boumaud 
et al. 2011), may drive energetic outflows expelling gas out 
of the galaxy and heating the halo. In addition, a so-called 
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“radio-mode” feedback may suppress the cooling of gas onto 
the galaxy over a longer timespan (Croton et al. 2006). 

Moreover, morphological quenching, proposed by Martig 
et al. (2009) may switch off or reduce the efficiency of star 
formation in galaxies through the existence of a dominant bulge 
that stabilizes the gas disk against gravitational instabilities (see 
Saintonge et al. 2012 and Crocker et al. 2012 for observational 
hints in the local universe, and Genzel et al. 2014 for galaxies 
at high redshift). Another process proposed to be responsible 
for the shut down of galaxies is halo mass quenching (e.g., 
Birnboim & Dekel 2003; Keres et al. 2005). In this scenario, 
dark matter haloes exceeding a critical mass of M ~ 10 12 M Q 
are able to stop the flow of incoming cold gas onto their central 
galaxies via virial shock heating, which leads to a decrease 
of star formation and/or suppresses star formation over longer 
timescales. However, at 7 > 7cm (with 1 < 7 cnt < 3), gas is 
predicted to penetrate to the central galaxy through cold streams, 
even in massive halos (e.g., Dekel et al. 2009). 

The above quenching mechanisms could plausibly leave 
an imprint on the structure of galaxies, because they may 
lead to high central concentrations due to internal processes 
and/or major mergers. This idea triggered a plethora of studies 
investigating the relation between star formation activity (or 
absence thereof) and galaxy structure. Local galaxy surveys 
demonstrated that quiescence is strongly linked with structural 
and morphological parameters, such as the Sersic index, stellar 
mass density, and the central velocity dispersion (Kauffmann 
et al. 2003, 2006; Schiminovich et al. 2007; Bell 2008; Fang 
et al. 2013; Cheung et al. 2012). This suggests that a high stellar 
mass surface density in the center of a galaxy is connected to its 
quenching. 

Recently, high-resolution imaging facilities on board the 
Hubble Space Telescope ( HST ), including images from the NIR 
Wide-Field Camera 3 (WFC3), have allowed us to trace the 
relation between structure and stellar populations further back 
in cosmic time. Early studies on the basis of high-resolution 
HST imaging found a correlation between color and bulge 
fraction for large samples of galaxies at 7 ~ 0.7-1 (Bell 
et al. 2004; Weiner et al. 2005; Koo et al. 2005). The Cosmic 
Assembly Near-infrared Deep Extragalactic Legacy Survey 
(CANDELS; Grogin et al. 2011; Koekemoer et al. 2011) has 
played a particularly pivotal role, due to its unique combination 
of multi-wavelength high-resolution imaging, large sample size, 
and depth. A common theme of several of the early CANDELS 
studies is that a correlation between galaxy structure and stellar 
populations (a “Hubble sequence”) has already been in place 
since at least z ~ 2.5 (Wuyts et al. 2011; Bell et al. 2012; 
Wang et al. 2012; Lee et al. 2013; see also Cheung et al. 2012). 
Additionally, Bruce et al. (2012) show that the majority of SFGs 
have disk-dominated light profiles, and passive galaxies have 
bulge-dominated profiles directly from decomposed bulge+disk 
light-profile fitting of M > 10 11 M Q galaxies at 1 < z < 3. 
However, it should be noted that they also find a significant 
fraction (~30%) of passive disk-dominated galaxies and star- 
forming, bulge-dominated systems in this high-mass regime. 
These results have set on a firmer footing the pioneering studies 
by Franx et al. (2008) and Kriek et al. (2009), that used ground- 
based data and significantly smaller samples, respectively. At all 
epochs beginning with at least 7 ~ 2.5, SFGs follow a relatively 
tight relation of 0.2-0. 3 dex scatter between their ongoing star 
formation rate (SFR) and assembled stellar mass, with a zero 
point that increases with look-back time. This relation has been 
dubbed the “main sequence” (MS) of SFGs (Noeske et al. 


2007; Elbaz et al. 2007; Daddi et al. 2007; Rodighiero et al. 
2011). Typical SFGs that lie on the MS are best approximated 
by exponential disk profiles, while the quiescent galaxy (QG) 
population below the MS is better described with higher Sersic 
indices (Wuyts et al. 2011; Bell et al. 2012). Common to all 
of the above studies, however, is that they were carried out 
on monochromatic (albeit mostly rest-optical) imaging, often 
exploiting only a subset of the entire CANDELS data set. If 
a galaxy component (e.g., disk) forms significantly more stars 
than another component (e.g., bulge), its relative weight in light 
can differ significantly from its contributions in more relevant 
physical units of stellar mass. 

Deep panchromatic high-resolution imaging data sets enable 
us to go beyond measurements of light profiles, by reconstruct- 
ing the distribution of stellar mass (see, e.g., Zibetti et al. 2009 
for a detailed application to nearby galaxies). While pioneer- 
ing work on pre-existing high-7 data sets using the NICMOS 
camera on board HST was limited to small sample sizes (e.g., 
Elmegreen et al. 2009; Forster Schreiber et al. 2011), Wuyts et al. 
(2012) made use of resolved seven-band stellar-population syn- 
thesis modeling on a pixel-by-pixel basis of over 600 massive 
( M > 10 10 M e ) SFGs at 0.5 < 7 < 2.5 and found that galaxies 
are overall smoother and more centrally concentrated in mass 
than they appear in light. Caused by a combination of internal 
extinction and star formation history (SFH) variations (see also 
Guo et al. 2012; Lanyon-Foster et al. 2012; Szomoru et al. 2013; 
Wuyts et al. 2013), the presence of mass-to-light ratio variations 
within galaxies emphasizes the importance of making measure- 
ments on mass maps when characterizing the stellar structure 
of high-redshift galaxies. Guo et al. (2011) furthermore demon- 
strated that QGs at high redshift also feature color (and hence 
mass-to-light ratio) gradients. 

In order to further shed light on the connection between 
galaxy structure and quenching, we reconstruct stellar mass 
maps of a large mass-selected sample (> 10 10 M 0 ) of galaxies 
at 0.5 < 7 < 2.5 and subject those to a detailed structural 
analysis. For this purpose, we exploit the available multi- 
wavelength data in all five CANDELS /3D-HST fields with 
accurate redshifts from the 3D-HST survey (Brammer et al. 
2012; van Dokkum et al. 2011). We do not restrict the structural 
measurements to one-component (Sersic) fits, but rather explore 
two-component (bulge + disk) decompositions. The latter allows 
us to carry out a more direct comparison to semi-analytic models 
(SAMs; specifically those by Somerville et al. 2008, 2012, 
further developed by L. A. Porter et al. 2014, in preparation), 
whose prescriptions are formulated in units of disk and bulge 
components. 

The paper is structured as follows. In Section 2, we give 
an overview of the observations and sample. The construction 
of stellar mass maps through resolved stellar population mod- 
eling, and the methodology to derive structural parameters, is 
described in Section 3. We present our observational results in 
Section 4, followed by a comparison with SAMs in Section 5. 
Finally, we summarize and discuss the implications of our find- 
ings in Section 6. 

Throughout this paper, we quote magnitudes in the AB 
system, assume a Chabrier (2003) initial mass function 
(IMF), and adopt the cosmological parameters (<T W , il,\, h) = 
(0.3, 0.7, 0.7). 

2. DATA AND SAMPLE SELECTION 

The core data set used for this work is the deep space-based 
HST imaging from the CANDELS multi-cycle treasury program 
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(Grogin et al. 2011; Koekemoer et al. 2011), complemented 
with redshift information from the 3D-HST grism survey (van 
Dokkum et al. 2011; Brammer et al. 2012). 14 

In order to determine galaxy-integrated masses and SFRs, 
we make use of additional multi-wavelength data in the 
CANDELS/3D-HST fields, including space-based photometry 
from Spitzer/ 1 R AC, Spitzer/MIPS, and llerschel /PACS , and an 
array of ground-based facilities (see Skelton et al. 2014 for a de- 
tailed description of the //-selected multi-wavelength catalogs 
produced by the 3D-HST team). 

The galaxy-integrated properties are derived following iden- 
tical procedures as Wuyts et al. (2011). These include stellar 
masses, based on U- to - 8 fim broadband spectral energy distri- 
bution (SED) modeling using population synthesis models from 
Bruzual & Chariot (2003) and SFRs derived from a “ladder of 
SFR indicators.” The latter method uses detected emission in ei- 
ther UV + PACS for PACS-detected galaxies (Lutz et al. 2011; 
Magnelli et al. 2013) or UV + MIPS 24 /im for MIPS-detected 
galaxies to compute the sum of the obscured and unobscured 
SFR. For galaxies lacking an infrared detection, the SFR is 
adopted from the best-fit SED model. We verified that consis- 
tent results are obtained when splitting star-forming and QGs 
based on the above measures of star formation activity or when 
adopting the rest-C/V7 diagram (Wuyts et al. 2007; Williams 
et al. 2009) 

2.1. HST Imaging 

The HST CANDELS observations used for this study 
comprise high-resolution imaging in five distinct fields: 
GOODS-South, GOODS-North, COSMOS, UDS, and EGS, 
covering a total area of 625 arcmin 2 . Typical limiting depths 
in //i6o are 27.0 mag for CANDELS/wide and 27.7 mag for 
CANDELS/deep (the central halves of the GOODS fields), 
respectively. Additional data used for this work include pre- 
existing Advanced Camera for Surveys imaging in the GOODS, 
EGS, and COSMOS fields (Giavalisco et al. 2004; Davis et al. 
2007; Koekemoer et al. 2007), plus WFC3 imaging in the Early 
Science Release (ERS; Windhorst et al. 2011). The available 
passbands are B 435 , V 60( , i 115 , z 8 so, his- H m for the GOODS 
fields, and V 606 , /si 4 , / 125 , // 16 O for the remaining fields. Addi- 
tionally, To 98 imaging was available for ERS as part of GOODS- 
South and T 105 for the regions with CANDELS /Deep coverage. 

All imaging used in our analysis was drizzled to a O'.' 06 pixel 
scale. For details on the observations and data reduction, we 
refer the reader to Koekemoer et al. (2011) and Grogin et al. 
( 2011 ). 

2.2. Sample Definition 

The main focus of this paper is on the structural shape of 
stellar mass distributions in a mass-selected sample of galaxies 
in the redshift range 0.5 < z < 2.5. To this end, we apply a 
mass cut M > 10 10 M G to select our galaxies, which is well 
above the mass completeness limit of the //^o-selected catalogs 
for all five fields. Our sample consequently spans a wide range 
of SFRs — from normal SFGs on the MS to QGs below — that 
already formed the bulk of their stars. For selection, we use the 
SED modeled galaxy-integrated mass estimates, but note that 
they nicely match the masses obtained by summing the resolved 
stellar mass distributions (see also Wuyts et al. 2012). 


14 The AGHAST G141 grism observations within GOODS-North from the 
GO-1 1600 program (PI: B. Weiner) are included in the analysis under 
“3D-HST” (Brammer et al. 2012). 


In order to determine redshifts for the galaxies in our sample, 
we use ground-based spectroscopic redshift information when- 
ever available. Otherwise, redshifts are fitted to the combination 
of 3D-HST grism data and broadband photometry. 

The total sample comprises 6764 galaxies, of which 3839 
and 2925 lie within the redshift range 0.5 < z < 1.5 and 
1 .5 < z < 2.5, respectively. In the following, these two redshift 
ranges are referred to as z ~ 1 and z ~ 2, respectively. 

3. METHODOLOGY 
3.1. Resolved SED Modeling 

A detailed description of the resolved SED modeling can be 
found in Wuyts et al. (2012). Here, we review only the key steps 
involved and additional processing steps. 

First, all images in the available wavelength bands are 
brought to the WFC3 Z /160 resolution (O'.' 1 8) by using the IRAF 
PSFMATCH algorithm. Next, a Voronoi-binning scheme from 
Cappellari & Copin (2003) is applied in order to ensure a 
minimum signal-to-noise ratio level of 10 for each bin in the 
corresponding H- band image. The binned multi-wavelength 
images are then fit with stellar population synthesis models 
from Bruzual & Chariot (2003), assuming a Chabrier (2003) 
IMF, a uniform solar metallicity, and a Calzetti et al. (2000) 
reddening law with visual extinctions in the range 0 < A v < 4. 
The adopted SFHs are exponentially declining, allowing for 
e-folding timescales down to 300 Myr. This choice is simplistic 
and may not be representative of the real SFHs (e.g., for 
SFGs, see Maraston et al. 2010). However, we are interested 
in the stellar mass, which, among SED-derived properties, is 
most robust against variations (and uncertainties) in model 
assumptions (e.g., Papovich et al. 2001; Shapley et al. 2005; 
Maraston et al. 2010). 

In order to conduct further structural analysis using the 
(binned) output of this SED modeling technique, we constructed 
pixelized (instead of Voronoi binned) images of the stellar mass 
surface density distribution. First, we construct a M/L map 
that, within the galaxy’s H\(,q segmentation map, has uniform 
values for pixels belonging to the same Voronoi bins. Pixels 
outside the segmentation map (i.e., containing sky noise and 
possibly faint wings of the galaxy extending below the signal- 
to-noise threshold) are assigned the average M/L of the nearest 
three Voronoi bins. The resulting expanded M/L map is then 
combined with the // 1 60 image to construct a final mass map 
at full //160 resolution. This ensures a smooth transition of 
the galaxy’s mass profile from the brighter central parts to the 
faintest regions. 

3.2. Structural Parameters 

In order to conduct a structural analysis of our galaxy sample, 
we employ the GALFIT (Peng et al. 2010) morphology fitting 
code to fit two-dimensional (2D) parametric models to the stellar 
mass distribution. For comparison purposes, we model the 2D 
//160 surface brightness distributions as well, following identical 
procedures. 

As part of our set of models, we use single Sersic models to 
parameterize the galaxy’s shape with the effective radius ( R e ) 
and Sersic index n. R e is the galactocentric distance containing 
half the total light/mass, and n is a measure of the cuspiness 
of the overall light/mass profile. We allow n to vary within the 
range 0.2 < n < 8. 

While the Sersic index n is often used and indeed serves as 
an approximate measure of the contribution of the bulge, it is 
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important to note that n does not translate one to one to the bulge- 
to-total ratio, B/T (see Appendix A; also Andredakis et al. 
1995; de Jong et al. 2004; Cibinel et al. 2013; Bruce et al. 2012; 
Brace et al. 2014). For a given bulge + disk composite profile, 
the best-fit Sersic index can be increased both by increasing 
the B/T , and by leaving the B/T constant while growing the 
extent of the disk relative to that of the bulge (i.e., lowering 
Re,B/ R e ,D)- For this reason, and to facilitate a more direct 
comparison to (semi-analytic) models that are parameterized 
in terms of bulge and disk units, we furthermore decompose 
the galaxies by fitting two-component (bulge + disk) models. In 
addition, we also performed one-component pure disk and pure 
bulge fits by forcing n = 1 and n — 4, respectively. 

For the bulge-to-disk decomposition, we adopt a procedure 
similar to that implemented in Brace et al. (2012) for the 
bulge+disk decomposition of //i6o light profiles, where we fixed 
the Sersic indices to be n — 1 for the disk and n — 4 for the 
bulge. 13 The centers are left free, but we restrict the relative 
distance between the bulge and disk centers to be less than 
2 pixels. All other parameters defining the two components 
( R e , the axial ratio b/a, and the total magnitude/mass of both 
components) are allowed to vary independently. 

In two-component modeling, the higher number of degrees of 
freedom increases the odds of the fit being trapped in a local / 2 
minimum. In order to mitigate this risk, we ran GALFIT using a 
grid of initial starting values. Our grid was constructed by using a 
range of size ratios between bulge and disk (R e ,B / Re.D), ranging 
from 0.1 to 1, in steps of 0.1. For each initial guess of R e b / R e ,D , 
the corresponding initial guess on B/T was then set, such that 
the Sersic index matching this initial configuration matches the 
one measured in the single-component fit (see Appendix A). 
Likewise, the initial magnitudes and absolute values of the 
initial size guesses for bulge and disk were set such that the 
total magnitude and half-light/mass radius of the composite 
profile matches the respective values determined from Sersic 
fitting. 

In cases where GALFIT yields solutions with implausibly 
small bulge sizes (<0.1 pixel, corresponding to <1/30 of the 
resolution, or <50 pc at z. ~ 2) and flags the outcome as 
potentially not converged and unphysical, we excluded that run 
from the grid. Also, solutions yielding a disk smaller in size 
than the bulge (i.e., R e .B/R e ,D > 1) were not included. 16 

After performing the fits for each point of the grid with 
initial guesses, we assigned a final B/T ratio for each object 
as the solution of the fit with the lowest x r e d - Here, the pure 
disk and pure bulge fits were also included. Their solutions 
generally show higher x r e d than the two-component fits, but are 
occasionally preferred over those in a x^j sense (oftentimes, 
these are bulgeless systems with n < 1). The two-component 
decompositions are statistically preferred over the single Sersic 
models (as based on both their x r ed va lues and the Akaike 


15 Since the choice of a Sersic index value for the bulge profile has long been 
debated in the literature (see Kormendy & Kennicutt 2004 and references 
therein), we repeated our two-component fits, with the only change fixing 

n = 2 for the bulge component to address the impact on the derived B/T 
values. We find a median difference of B n -A / T — B n =i/T = 0.03*/ / for 
QGs and 0.01//® for SFGs, with the errors marking the lcr scatter. Any 
systematic trends are small, and our results are therefore robust against the 
precise value of «bulge adopted. 

16 We note that nuclear stellar disks in early-type galaxies do exist (e.g., Jaffe 
et al. 1994; van den Bosch et al. 1994; Ferrarese et al. 1994), but they are 
impossible to resolve at the HST resolution for z ~ 2 galaxies. 


information criterion, AIC 17 ) for ~2/3 of the total sample. 
Those systems for which the single Sersic model yields a lower 
value of x r ed an( J AIC typically feature shallow profiles with 

n < 1. 

In both single- and two-component fitting, we use an au- 
tomated scheme that pre-determines neighboring sources that 
need to be masked or fitted simultaneously, and passes initial 
guesses of fitting parameters to GALFIT. Those include esti- 
mates on size, total magnitude, and center. The initial guess on 
size is based on the distance from the center to the radius at 
which the curve of growth reaches 50% of the galaxy’s total 
flux. The mass-weighted center of the galaxy derived within its 
segmentation map is adopted as initial estimate for the center. 
GALFIT takes into account the convolution of the model with 
the point-spread function (PSF). Both for fitting light and mass, 
we use a PSF that is a combination of stacked stars and a Tiny- 
Tim (Krist 1995) model PSF. For a more detailed description of 
the used PSF, see van der Wei et al. (2012). 

We emphasize that throughout this paper our working def- 
inition adopted for bulge and disk components is based on 
the above bulge+disk modeling of stellar mass or light maps, 
as there are currently no empirical constraints on whether or 
not stars assigned to a bulge/disk component are dynamically 
hot/cold. As the SINS survey of z ~ 2 galaxy kinematics 
demonstrated that high-redshift SFGs are dynamically (Forster 
Schreiber et al. 2009) and morphologically (Forster Schreiber 
et al. 2011) distinct from local spiral galaxies, we caution that 
bulge fractions as derived by our decompositions may have a 
somewhat different meaning than they would have in the local 
universe. 

In addition, we estimated a typical measurement uncertainty 
on B/T by setting up an array of model galaxies with a range 
of mass, size, B/T ratio, and R e ,B/Re.D, which is similar to 
our data. The grid of model galaxies consisted of 5, 3, 11, and 
10 grid points in mass, R e , B/T, and R e ,s/Re,D , respectively. 
These were then inserted in multiple empty sky regions of the 
CANDELS UDS field to mimic the typical background noise. 18 
We next ran GALFIT on all 8250 mock galaxies using our two- 
component fitting scheme. The measurement error on B/T of 
each galaxy in our sample given its mass, radius, and profile 
shape is then finally assigned as the scatter among the recovered 
B/T ratios of the corresponding model galaxy. 

Typical measurement errors in B/T for star-forming and QGs 
are, on average, ~0.05 and ~0.06 at z ~ 1 and ~0. 1 and ~0. 13 
at z ~ 2, respectively. The distribution of errors peaks below 
the median (<0.05), and shows a tail toward higher B/T errors. 
Two alternative methods to estimate the uncertainty in B/T, 
namely, the formal random uncertainty reported by GALFIT 
and re-fitting the observed galaxies after applying additional 
background and Poisson noise, generally lead to lower estimated 
uncertainties (by a factor of ~2 in the case of GALFIT). In the 
remainder of the paper, we therefore adopt the most conservative 
error estimates inferred from our analysis of the inserted mock 
galaxies. 


17 In evaluating models with the AIC, the preferred model is the one that 
minimizes / 2 + 2 p + 2 p(p + 1 )/N — p — 1), where p is the number of free 
parameters in the fit and N is the number of data points used in the fit. 

18 UDS is part of CANDELS-Wide, which was exposed for two HST orbits 
divided over F125W and F160W. Part of our data set comes from the 
CANDELS-Deep regions (the central halves of the GOODS fields), which 
received four orbits per pointing in F125W and F160W each. The inferred 
uncertainties from our analysis of mock galaxies can therefore be considered 
as conservative estimates. 
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4. RESULTS ON GALAXY STRUCTURE 
4.1. The Evolving Mass Budget of Disks and Bulges 

When exploiting the bulge-disk decompositions of the stellar 
mass maps derived for our sample of massive galaxies, we first 
evaluate the average mass budget of disks and bulges. The 
following is an example of picking a random star out of the 
sample of massive galaxies above 10 10 M Q . At 1.5 < z < 2.5, 
the probability that this star belongs to a bulge component is 
46%. Increasing the mass limit to log(M) = 10.5 or 11 yields 
a higher probability for the star to be associated to the bulge 
(49% and 54%, respectively). Perhaps somewhat surprisingly, 
the fraction of stars residing in a bulge component rises only 
slightly, to 0.5 < z < 1.5 (47%, 50%, and 56%) for galaxies 
more massive than log(M) = 10, 10.5, and 1 1, respectively. 

The formal uncertainties to the above stated probabilities, 
including sample variance and typical measurement errors on 
B/T, are limited to a few percent. The total error budget 
is likely dominated by systematics, for example, related to 
the assumptions made in stellar population modeling (see 
Section 3.1). We note, however, that only M/L uncertainties 
with a differential impact on bulges and disks will affect 
the above numbers. Even if the M/L ratio of bulges were 
systematically under- or overestimated by 0.2 (0.3) dex with 
respect to those of disks, the change would be limited to 
~7(10)%. 

Since bulges, unlike stellar disks, can be considered sinks in 
the continuous assembly of a galaxy’s stellar component, the 
rising mass density of stars in bulges (by a factor of ~1 .8 from 
the higher to the lower redshift bin) seems to be compensated 
largely by the continuing assembly of new disks. Splitting our 
sample in liner redshift bins, we do find the fraction of stellar 
mass in bulges to increase more significantly, by a factor ~1.5 
over the entire 6 Gyr timespan sampled in our study. 

Overall, the bulge mass fractions are higher than what would 
be inferred from fits to the //-band surface brightness profiles, 
as the median mass-to-light ratio of disk components is 0.2 dex 
lower than the bulge components. The above numbers address 
the evolving mass budget of disks and bulges for a mass-limited 
sample, including both star-forming and QGs. In the remainder 
of the paper, we will break down our sample by star formation 
activity. 

4.2. Profile Shape 

In recent years, several //ST-based studies have investigated 
the structural differences between star-forming and QGs at high 
redshift. In all of these analyses star-forming systems are found 
to have significantly larger rest-optical sizes than their quiescent 
counterparts at the same mass and redshift (e.g., Toft et al. 2009; 
van der Wei et al. 2014). In addition, their surface brightness 
profile shapes tend to be shallow (n ~ 1), whereas QGs feature 
cuspier light profiles (e.g., Wuyts et al. 2011; Bell et al. 2012; 
Cheung et al. 2012, although see also Bruce et al. 2012 for 
decomposed light profiles at high stellar masses and redshifts). 

At the same time, the same multi-wavelength, high-resolution 
look-back surveys have also established that substantial mass- 
to-light ratio variations in the rest-optical can occur, not only 
between but also within galaxies (see, e.g., Wuyts et al. 2012; 
Guo et al. 2012; S. Boada et al. in preparation). Since new stars 
typically tend to form from gas settled in a disk configuration, 
M/L ratio variations are anticipated to give rise to a composite 
light profile in which the disk component has a relatively larger 
weight (per unit mass) than the bulge. It is therefore important to 


address to what degree the above structural distinction between 
the two classes of galaxies is intrinsic to their distribution 
of stellar mass, or, conversely, can be attributed to stellar 
population effects. The answer is of immediate relevance to 
our understanding of quenching, as the latter scenario could 
imply that compact quiescent systems can evolve from the star- 
forming MS by simple fading, without invoking an associated 
morphological transition. Kriek et al. (2009) investigate this 
scenario for a spectroscopically confirmed sample of massive 
Z ~ 2 galaxies, finding that three out of six massive star-forming 
systems have dense cores, and thus may passively evolve into 
compact galaxies due to the fading of the outer star-forming 
regions. Szomoru et al. (2010), on the other hand, exploit the 
exquisite depth of the Hubble Ultra Deep Field to probe the 
surface brightness profile of a massive compact QG at z = 1.91, 
ruling out the existence of a faint extended envelope or disk 
around the observed galaxy. Another argument against fading 
comes from Cheung et al. (2012), who derived the stellar masses 
of bulges for both star-forming and QGs at 0.5 ^ z < 0.8. They 
found that the bulges of SFGs are half as massive as those of 
similar-mass QGs, implying they cannot simply fade onto the 
red sequence without structural evolution. 

Using the stellar mass maps reconstructed for our mass- 
selected sample of 6764 galaxies at 0.5 < z < 2.5 with 
log (M) > 10, we are now able to draw statistically significant 
conclusions on the structural distinction between high-z galaxies 
prior to and after quenching. In Figure 1, we compare the 
shape of the stellar mass distributions (i.e., corrected for spatial 
M/L variations) of star-forming and QGs, and study their 
dependence on the total galaxy stellar mass. We consider profile 
parameters based on single-component (i.e., Sersic) fits as well 
as two-component (bulge + disk) decompositions, and show the 
results for two separate redshift intervals: 0.5 < z < 1.5 and 
1.5 < z < 2.5. In both cases, we identified a galaxy as quiescent 
if its specific SFR (sSFR) satisfied sSFR < l/3fHubbie, where 
htubbie is the Hubble time at the redshift of the galaxy, and as 
star-forming otherwise. We tested that a definition of quiescence 
based on the location of a galaxy in the UVJ diagnostic diagram 
(Wuyts et al. 2007; Williams et al. 2009) yields effectively 
identical results. 

Figure 1 immediately highlights that the distinct structural 
appearance of star-forming and QGs is intrinsic to its internal 
distribution of stellar mass, and not just driven by stellar 
population or obscuration effects. In fact, a comparison to the 
equivalent plots based on //-band surface brightness profiles 
rather than mass maps (see Appendix B) indicates that stellar 
population effects (when measuring at rest-optical wavelengths) 
only induce a modest, albeit non-negligible shift. At all masses, 
QGs feature cuspier stellar mass distributions (i.e., higher n ) 
than star-forming systems. Their typical best-fit Sersic index is 
also an increasing function of galaxy mass. Interestingly, among 
SFGs the profile shape is not independent of stellar mass. An 
increase in n is apparent above 10 11 M Q , both at z ~ 1 and at 
Z ~ 2. A similar trend of increasing cuspiness at the tip of the 
MS was noted by Wuyts et al. (2011, Figure 1; see also E. J. 
Nelson et al. in preparation). 

Next, it is worthwhile to reflect on what is measured when 
fitting Sersic profiles. Appendix A illustrates that, when consid- 
ering galaxies as superpositions of bulge and disk components, 
a given best-fit Sersic index does not necessarily correspond 
one-to-one to a unique B/T value, even though it is often inter- 
preted as such. Given a bulge+disk system with an associated 
best-fit n, one can increase its n either by boosting B/T , or. 
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Figure 1 . Sersic index and B /T ratio as a function of the total stellar mass of our galaxy sample spanning the redshift range 0.5 < z < 2.5. The whole sample (black 
non-filled symbols) is divided in star-forming (blue symbols) and quiescent galaxies (red symbols). The respective shaded areas mark the 50th percentile scatter of 
the distributions within one bin, while the error bars indicate the uncertainty on the median value. Results for the lower and upper halves of the redshift interval are 
marked with dotted and dashed lines, respectively. The bins contain (from left to right) 1439, 1091, 886, 693, 484, 259, 1 14, and 35 star-forming systems, and 227, 
267, 298, 369, 280, 214, 83, and 25 quiescent systems within the whole redshift range. SFGs show a clear trend of increasing Sersic index and B/T with increasing 
stellar mass, reaching B/T ~ 0.4-0. 5 above 10 11 Mq and reflecting the build-up of central mass concentrations in main-sequence SFG up to z ~ 2.5. 

(A color version of this figure is available in the online journal.) 


alternatively, by growing the extent of the disk with respect to 
that of the bulge without any change to B/T. 

The right-hand panel of Figure 1 explores the B/T ratio as a 
function of galaxy mass, for SFGs and QGs separately. Again, 
we find a clear anti-correlation between star-formation activity 
and bulge prominence. Focusing on the star-forming population, 
the median B/T is limited to below 30% for intermediate 
mass SFGs (10 < log(M) < 11), while typical bulge mass 
fractions rise to 40%-50% above 10 11 M Q . We note that there 
is a significant scatter in the distribution of individual B/T 
values around the median for both quiescent and SFGs. We 
investigated the variation in median trends when varying the 
binning intervals, finding negligible changes at lower masses, 
while the median B/T of the most massive (log(M) ~ 11.3) 
SFG bin changes by ±0.1, depending on the applied binning 
intervals. 19 

We note that measurements on the H band yield bulge mass 
fractions among SFGs that are lower by ~30% (on average), 
as can be understood from a disk component composed of a 
younger, lower M/L stellar population than the bulge. 

From the two-component fits, we infer a typical R e ,B/Re,D 
size ratio of ~0.2, albeit with significant scatter (see Figure 6). 
The median size ratio shows little dependence on star-formation 
activity or mass, over the range probed by our sample. Given 
the enhanced B/T values in QGs, and the fact that bulges 
have smaller half-mass radii than disks, one might wonder if 
the difference in total size between SFGs and QGs can be 
completely accounted for by a redistribution of stellar material 
from the disk to the bulge, without changing the extent of each 
of the components individually. Our analysis confirms that the 


19 The binning scheme applied in Figure 1 is such that the most massive bin 
still contains more than 10 galaxies, allowing a robust estimation of the 
median. 


change in the B/T of SFGs prior to or during quenching is, 
to a large extent, responsible for the size difference between 
the quiescent and star-forming population. However, some 
fraction of the shrinking size is still attributed to the individual 
components being smaller. In Figure 2, we show the total sizes as 
well as the sizes of the individual components for star-forming 
and QGs, as measured on the mass maps. While the total sizes 
of SFGs and QGs are noticeably different (by a factor ~3 at 
log(M„) ~ 10.5), the difference in the size of the bulge and 
the disk components between SFGs and QGs, respectively, is 
smaller, typically by a factor ~1.5. 

4.3. Fraction of Quenched Galaxies 

With the morphological parameters of the mass maps for our 
entire galaxy sample in hand, we now proceed to relate those 
with galaxy-integrated star formation properties. 

The three panels of Figure 3 show, from left to right, the 
fraction / que nc h of quenched galaxies as a function of total 
stellar mass, bulge mass, and disk mass, respectively. Here, 
we again define galaxies as quenched/quiescent when sSFR < 
1/3 * tHubbie(t(), and as star-forming otherwise. 

The uncertainties in ./quench are derived via a bootstrapping 
method and represent the 68% confidence levels. They include 
sample variance and the typical measurement errors on B/T, 
which are derived as described in Section 3.2. For the boot- 
strapping, we computed / quen ch for 1000 samples, which are 
randomly drawn from the original sample, with replacement. 
For each bootstrap iteration, we displace the B/T values for 
each galaxy by the typical measurement error in B/T, given the 
galaxy’s magnitude, size, and measured profile shape. 

The top left panel of Figure 3 illustrates that the fraction of 
quenched galaxies increases with increasing mass, from ~0. 1 at 
around 10 10 M 0 to ~0.5 at 2 x 10 11 M Q . A second conclusion 
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Figure 2. Effective total size, bulge size, and disk size as a function of total stellar mass in the redshift range 0.5 < z < 2.5. The total sizes are derived from single 
Sersic fits, and bulge and disk sizes are inferred from the bulge-disk decompositions. The whole sample (black symbols) is divided into star-forming (blue symbols) 
and quiescent galaxies (red symbols). The respective shade areas mark the 50th percentile scatter of the distributions, while the error bars indicate the uncertainty on 
the median value. There is a clear difference in total size between SFGs and QGs over the whole mass range. In addition, the bulge and disk components individually 
appear to be smaller for the star-forming population up to masses of log (M) ~ 1 1.3. 

(A color version of this figure is available in the online journal.) 


is that the fraction of quenched galaxies is overall higher for 
the lower redshift bin, by, on average, a factor of ~2. Both of 
these results are well established in the literature. The rising 
mass function of QGs over cosmic time has most recently been 
quantified on a firm statistical footing by Muzzin et al. (2013) 
and Ilbert et al. (2013), both of which exploit the wide-area 
UltraVISTA survey. What CANDELS lacks in number statistics 
compared to UltraVISTA, it adds in depth and high resolution. 
Exploiting these key strengths, we now turn to the dependence 
of the quenched fraction on galaxy subcomponents: the mass of 
their bulge (middle panel) and disk (right-hand panel). Clearly, 
the dependence of / que nch on the bulge mass is much stronger 
than on the disk mass, which does not show any significant 
correlation with / quenc h above log(MDisk) ~ 9.5. Toward lower 
disk masses, / que nch increases rapidly, but we point out that this 
trend is entirely driven by the (total) stellar mass limit of our 
sample (log(M*) > 10; that is, the galaxies occupying the lowest 
iff Disk bins are necessarily heavily bulge-dominated systems, 
which tend to form relatively few stars). If lower mass galaxies 
were to be included, less massive, disk-dominated SFGs would 
likely outnumber these massive spheroids with small residual 
disks in the low M D ^ bins, producing a flat relation of / que nch 
with Moisk over the full range probed. Above respective masses 
of 10 10 Mq, / q uench increases more rapidly with bulge mass than 
with total stellar mass in both redshift ranges (~0.35 per dex of 
Tf bulge compared to ~0.3 per dex M*, or ~0. 1 per dex of Moisk)- 
With the bulge-to-disk decompositions in hand, we next split 
the galaxy sample in bins of B / T , and explore second parameter 
dependencies. Considering first the dependence of / que nch on 
the total stellar mass, it is apparent that, at a given total mass, 
/quench is increasing significantly with increasing B/T ratio. 
The middle panels of Figure 3 illustrate that, when considering 
the dependence of / quen ch on bulge mass, the different B/T bins 
align along a much tighter locus. In contrast, a large spread is 
seen as a function of disk mass (right-hand panels of Figure 3). 
In order to quantify these trends, we compute the Spearman’s 
rank correlation coefficient ( r s ) for the relations of / quen ch 
with log(M„), log(MDisk), and log(Mb u i ge ) for respective masses 
log(M) >10. We find that r s is indeed significantly higher for 


the relation / que nch versus log(Mb u i ge ) (G ~ 0.68), than for both 
fq uench versus log (AT*) (G ~ 0.32) and f quench versus log(M Disk ) 
(r s ~ —0.05), as measured for z ~ 1. Consistent results are 
found for z ~ 2. 

We investigated the impact of defining quiescence based 
on a UVJ color-color criterion instead of a sSFR cut. When 
applying a CVT-based selection of QGs, we find an overall good 
agreement with the trends presented in Figure 3. Quantitatively, 
small changes occur, with / que nch increasing by ~7% for the 
entire 0.5 < z < 2.5 sample integrated over all masses. The 
good agreement is not surprising, since the precise threshold 
in sSFR used to select QGs (sSFR < 1/3 * /Hubbie(Q) was 
chosen to yield maximum overlap with the UVJ selection 
criterion. 

Taken together, this demonstrates that the build-up of a 
bulge seems to play a critical role in the quenching process 
of galaxies, whereas the disk does not. The amount of stars in 
the disk component of a galaxy has little to no predictive power 
regarding its star-forming or quenched state, unless B/T (and 
hence the bulge mass) is also known. We find a qualitatively 
similar behavior at z ~ 2 as at z ~ 1 , but note that the cosmic 
evolution in the quiescent fraction cannot solely be attributed to 
continuing bulge growth over time, as galaxies in the same M* 
and B/T bin at z ~ 1 are more likely to be quenched than those 
at z ~ 2. Appendix B illustrates how the equivalent diagrams 
composed from fits to the //-band surface brightness, rather than 
the stellar mass distribution, exhibit a larger spread from low to 
high B/T bins. This generic behavior can be understood from 
a physical picture, where the disk component has a relatively 
larger weight in light than in mass. 

Our work is in agreement with, and takes the next step beyond, 
previous reports that the inner stellar mass density is better 
related to the SFH than the total stellar mass (Franx et al. 2008; 
Bell et al. 2012), as inferred from rest-optical imaging of smaller 
samples of high-redshift galaxies (see Kauffmann et al. 2003 and 
Fang et al. 2013 for a local universe reference, and Cheung et al. 
2012 for intermediate redshifts z < 0.8) 

Importantly, the same behavior explored here over the redshift 
range 0.5 < z < 2.5 extends in a strikingly similar fashion all 
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Figure 3. Fraction of massive galaxies ( M > 10 10 Mq) that are quiescent (/quench), as a function of total stellar mass (left columns), bulge mass (middle columns), 
and disk mass (right columns) for z ~ 1 and z ~ 2. Galaxies with sSFR < 1/3 * /Hubble are assigned as quiescent, the others are assigned as star-forming. In the middle 
and bottom rows, we split the galaxy sample in bins of B/T for the two redshift bins. Uncertainties are derived from bootstrapping and include sample variance as 
well as measurement uncertainties. A positive trend of /quench is seen with total stellar mass, with the trend becoming stronger when correlating /quench with the bulge 
mass, whereas /quench shows no positive correlation with the mass of the disk component. At a given galaxy mass, /q uen ch is increasing significantly with increasing 
B / T ratio. The scatter in /quench among the different B/T bins is largely reduced when correlating against the bulge mass, implying that the mass within the bulge 
of the galaxy is correlating best with quiescence. This trend is qualitatively similar for both redshift ranges, with the overall quenched fractions being lower at z ~ 2 
than at z ~ 1 . 

(A color version of this figure is available in the online journal.) 


the way to the present day, as demonstrated by Bluck et al. 
(2014) who exploit the large number statistics of SDSS. 

5. COMPARISON WITH SAMS 

5.1. The Somerville Model 

SAMs have a rich history of trying to reproduce galaxy scaling 
relations and abundances, with the goal of guiding our interpre- 
tation of the observational results. Here, we focus specifically on 
the SAM developed by Somerville et al. (2008) and further up- 
dated by Somerville et al. (2012) and L. A. Porter et al. (2014, 
in preparation), which is rooted in the Bolshoi cosmological 
dark matter simulation (Klypin et al. 201 1). :() As is generic to 

20 Hereafter, we refer to this model as the Somerville et al. SAM. 


all SAMs, the model relies on simplified analytic prescriptions 
for the dynamical and astrophysical processes down from entire 
galaxy scales, rather than on kiloparsec to parsec scales (the res- 
olution below which state-of-the-art cosmological and zoom-in 
hydro-simulations resort to subgrid physics, respectively). This 
limitation, however, yields the enhanced flexibility of a rela- 
tively inexpensive runtime, allowing the straightforward gener- 
ation of statistically significant model galaxy populations, and 
the tuning of parameters to observational constraints, such as 
mass functions and scaling relations (only empirical constraints 
from the nearby universe were used in tuning the parameters of 
the model considered here). The fact that SAMs conceptually are 
formulated in units of bulge and disk components furthermore 
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Figure 4. Same as Figure 3 using galaxies from the Somerville et al. SAM. Galaxies with sSFR > 1/3 * ^Hubble are assigned as star-forming, the others are assigned 
as quiescent. Uncertainties are derived via bootstrapping to reflect the sample variance. 

(A color version of this figure is available in the online journal.) 


makes them suitable for a direct and meaningful comparison to 
the diagnostics explored in this paper. 

A detailed description of the prescriptions for cooling, star 
formation, feedback, and structural growth is provided by 
Somerville et al. (2008, 2012), with extensions and applications 
to the CANDELS data set presented by L. A. Porter et al. (2014, 
in preparation). The input to and output from the model is further 
contrasted to that of other SAMs by Lu et al. (2013). 

For a detailed discussion of the physical recipes of this 
SAM and the resulting output in the context of a larger set 
of SAMs, we refer the reader to Lu et al. (2013). For the 
sake of the comparison presented here, we emphasize that 
none of the relations investigated in this paper were part 
of the set of observational constraints to which the model 
parameters were tuned. Model parameters of the SAM were 
tuned to (approximately) match the global stellar mass function, 
the stellar mass function of early- and late-type galaxies, the 
gas fraction as a function of stellar mass for disks, and the 
mass-metallicity relation for stars. 

Also of particular relevance is the fact that, in the model, 
stars form either in the disk following a Kennicutt-Schmidt 
law (Kennicutt 1998), where the disk scalelength is set by a 
similar methodology as Mo et al. (1998), or, in the event of a 
merger during a starburst. Bulge formation, as well as feeding 
of the central supermassive black hole, can happen through 
two channels: mergers or disk instabilities (see L. A. Porter 
et al. 2014, in preparation). The starburst, black hole accretion, 
and morphological transformation induced by mergers depends 
on the mass ratio and gas fraction, as calibrated using a large 
suite of binary merger simulations (Somerville et al. 2008, 2012; 


Hopkins et al. 2009). Star formation is moderated by heating the 
gas by supernovae as well as through AGN feedback. In addition 
to the quasar mode, during which AGN can drive powerful 
outflows, black holes also grow more gradually over longer 
timespans through the so-called radio mode (i.e., suppression of 
cooling via radio jets). No explicit connection between the bulge 
mass and quenching (as may for example be expected from 
the Toomre Q stability criterion in a gravitational quenching 
scenario; see Section 6.2) was built into the model. 

5.2. ^/quench U7 th(? SAM 

In Figure 4, we show the model equivalent of Figure 3, de- 
scribing how the fraction of galaxies that are quenched depends 
on the total stellar mass, the mass of stars in the bulge, and the 
mass of stars in the disk component, respectively. As for the ob- 
servations, we define the threshold for a galaxy to be quenched 
with a ruler moving with redshift: sSFR < l/(3fHubbie)- This def- 
inition, rather than an application of the UVJ diagnostic, remains 
closer to the direct output of the SAM, avoiding a translation 
to mock spectral energy distributions, which would introduce 
additional assumptions and uncertainties. 21 

At first glance, the SAM features several of the characteristic 
trends noted earlier for the CANDELS galaxies: / que nch rises 
toward later cosmic times, increases with the total stellar mass 
(more steeply so with bulge mass), and shows no appreciable 


21 It should be noted that equivalent assumptions and uncertainties associated 
with the conversion from light to physical properties enters upon SED 
modeling of the observed galaxies. The choice of how far to take the models to 
the observations, or vise versa, therefore remains somewhat arbitrary. 
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correlation with disk mass above log(M D i s i;) = 9.5. In the 
interval 9 < log(MDisk) 1$ 9.5, a sharp drop with increasing 
Moisk is noted, as also shown in the observations. Given 
the M > 10 10 M q threshold of our sample selection, the 
latter objects are necessarily heavily bulge-dominated. Without 
imposing such a mass limit, less massive, disk-dominated 
SFGs would outnumber massive early-type galaxies with small 
residual disks in the lower Moisk bins. Despite the qualitative 
success, quantitative differences in the quenched fractions of 
model galaxies with respect to those observed are clearly 
present. The discrepancy is most severe at z ~ 2, where modeled 
/quench values are, on average over the whole displayed mass 
range, of order a factor ~2-3 short of the observed, hinting at 
an underestimated quenching rate and/or inefficiency to prevent 
quiescent systems from rejuvenating. 22 

When splitting the SAM galaxies at each M„, Mb u i ge , and 
Moisk in bins of B/T (middle and bottom panels of Figure 4), 
we reproduce a similar behavior as found for the real universe 
in Section 4.3. Namely, the total stellar mass acts as a poorer 
predictor of the quenched state of a galaxy. This situation can 
be remedied if in addition to M*, B/T (and hence Mb u i ge ) is 
also known. Quantitatively, the correlation between / quen ch and 
Mbuige is measured to be the strongest ( r s ~ 0.46), whereas M* 
and Moisk only show weak correlation with / que nch ( r s ~ 0.21 
and r s ~ —0.1, respectively). The quoted values of r s are 
measured for respective masses of log (M) > 10 and at z ~ 1. 
The values for r s at z ~ 2 are similar, with the correlation 
between / que nch and log(M*), as well as between / que nch and 
log(Moisk), being somewhat stronger. At 0.5 < z < 1.5, less 
than 20% of all massive galaxies with B/T < 0.2 are classified 
as quiescent. Conversely, the majority of galaxies in the upper 
B/T bin (with B/T > 0.8) have low sSFR. These inferences 
are in agreement with the SAM and the observations. Also in 
agreement is the fact that Mb u i ge serves as a better predictor of 
/quench than the total stellar mass, with different B/T bins being 
more (albeit not perfectly) aligned along a single locus in the 
/quench versus M bu i g e diagram. 

At 1.5 < z < 2.5, the model predictions are skewed 
toward too low / que nch values, as noted earlier. However, in 
relative terms the same generic behavior, as a function of bulge, 
prominence is notable. 

We note that most of the trends in Figure 4 are driven by 
physical prescriptions in the SAM affecting central galaxies 
rather than satellites, as centrals account for 80% (90%) of 
the model galaxy population above log(M) = 10 (11). Those 
massive galaxies classified as satellites are further subjected to 
additional environmental quenching processes, resulting in a 
higher / que nch(M*, Mb u i ge ) for this particular subpopulation. 

5.3. The Agent of Quenching 

Given the qualitative agreement between the model and 
observations, we can now pose the question how, in the context 
of the Somerville SAM, the relation between structure and 
stellar populations could be interpreted. To this end, we consider 
the dependence of / que nch on two physical properties of galaxies 
in the SAM that are not observationally accessible for our 
CANDELS sample: the halo mass Mn a io and the mass of 
the central supermassive black hole Mbh- The top panels 
of Figure 5 illustrate that the probability of a galaxy being 
quenched increases toward high M Ha io and high Mbh- In detail, 


22 See also Ciambur et al. (2013) for a discussion on the quenched fraction in 
the Garching SAMs. 
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Figure 5. /quench as a function of halo and black hole mass using galaxies from 
the Somerville et al. SAM. In the middle and bottom rows, the galaxy sample 
has been split in bins of B/T . Galaxies with sSFR > 1/3 * ^Hubble are assigned 
as star-forming, the others are assigned as quiescent. Uncertainties are derived 
via bootstrapping to reflect the sample variance. 

(A color version of this figure is available in the online journal.) 


however, the dependencies on the two look different. The 
increase in / quen ch with Mnaio is gradual over nearly two orders of 
magnitude. In contrast, a much sharper upturn of / que nch emerges 
above log(MBH) = 7.5. This behavior is especially notable at 
Z ~ 1, but a rise above the same threshold is also present in the 
Z ~ 2 population. Breaking the model galaxy population down 
by its structural properties, we find a wide spread in / que nch 
for different B/T at a given /V/ Ha i 0 . A different, less scattered 
behavior would be expected if halo mass quenching were the 
sole and dominant mechanism (Birnboim et al. 2007; Dekel & 
Birnboim 2008; Dekel et al. 2009). However, as a function of 
Mbh, a similar upturn in / que nch is present for all B/T bins 
equally above log(MBH) = 7.5. 

The increased scatter in the relation between / que nch and 
Mb u i g e, compared to the tight correlation of / quenc h and Mbh can 
be explained by the SAM’s MBH-Mb u i ge relation. For a given 
B/T bin, the scatter in the MBH-Mb u i ge relation is significant 
compared to the dynamic range in the bulge masses plotted. 
At a given bulge mass, the scatter stems in part from an anti- 
correlation between the black hole mass and the level of star 
formation activity (sSFR). 

Our analysis illustrates that in the Somerville model, which 
includes feedback from both supernovae and AGN, the central 
supermassive black hole acts as the primary agent of quenching 
in massive galaxies, and its accumulated mass (i.e., the integral 
over past accretion activity) is tightly related to the probability of 
finding a galaxy in a quenched state. Since the physical processes 
giving rise to bulge and black hole growth are the same in the 
SAM (mergers and disk instabilities), the stronger relation of 
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/quench with Mb u i ge than with M*, present in both observations 
and model predictions, is not surprising. 

We stress that in the Somerville model, no direct causal link 
between the presence of a bulge and quenching is implemented. 
The bulge is simply the accessible observable that correlates 
most tightly with the actual agent of quenching in this par- 
ticular model: the supermassive black hole. Observationally, 
there is increasing evidence for AGN-driven outflows of mas- 
sive z ~ 1-3 galaxies (Forster Schreiber et al. 2014; Nesvadba 
et al. 2011; Harrison et al. 2012; Cano-Dfaz et al. 2012). Given 
the shortcomings of the SAM in a quantitative sense, notably 
its under-prediction of the quiescent population at z ~ 2, ad- 
ditional or other quenching processes may be at play in the 
real universe. One such process, which is causally linked to the 
presence of a bulge, could be morphological quenching (Martig 
et al. 2009; Genzel et al. 2014). In such a scenario, the high 
central stellar density provided by a bulge stabilizes the gas 
disk and prevents it from forming stars. While plausibly only a 
temporary measure (as no gas is expelled, or stopped from ac- 
creting through this mechanism), it could potentially contribute 
to suppressing star formation in z ~ 2 galaxies more efficiently 
and/or prevent them from returning to the star-forming branch 
in the SAM. In the local universe, star formation efficiencies of 
bulge-dominated systems are reduced by factors of ~2-3 com- 
pared to disk-dominated galaxies (Saintonge et al. 2012; Martig 
et al. 2013). 

6. DISCUSSION 
6.1. Structural Change 

In order to study the morphological differences between SFGs 
and QGs, and to draw conclusions on the possible structural 
changes of SFGs as they move along the MS, we first examined 
the mass dependence of the profile shape of SFGs and QGs as 
traced by the Sersic index and B /T ratio. We have shown that 
QGs are structurally distinct from the star-forming population as 
seen by overall higher Sersic indices and B /T ratios at a given 
stellar mass. SFGs show rising trends of their median Sersic 
index and B /T ratio with increasing stellar mass, with the latter 
rising up to ~40%-50% above 10 11 M 0 . These findings give 
insights about the link between the structural evolution of SFGs 
and the quenching process as they move along the MS. 

Analyzing the Schechter functional forms of the SFG and 
QG mass function as a function of redshift, Peng et al. (2010) 
conclude that the quenching rate of galaxies climbing the MS 
rises proportionally to the SFR (and given the near-linear 
MS slope therefore also proportionally to the stellar mass, 
hence their terminology “mass quenching”). This corresponds 
to a survival probability on the MS that drops exponentially 
with mass, implying that, while nearly all low-mass SFGs are 
destined to continue growing along the MS, toward the high- 
mass end the MS becomes progressively more dominated by 
nearly dead SFGs. In fact, the sub-unity slope of the MS, and 
possible flattening at the high-mass end (Whitaker et al. 2012), 
may well be interpreted in this context — the typical SFG above 
10 1 1 M 0 is already undergoing some level of quenching, thereby 
deviating from the projected path along a SFR-mass relation 
of slope unity, which could be expected from cosmological 
accretion rates in the absence of quenching. Tying in our 
observational results on galaxy structure, the deviation toward 
high median n and B / T at the massive end reflects the typical 
structure of nearly dead star-formers that account for the bulk of 
SFGs above 10 1 1 M 0 . The fact that they are structurally distinct 


implies that the morphological transition happens first, to be 
followed later by the departure from the MS. Bulge growth 
precedes quiescence. Such a morphological change prior to 
quenching is in line with qualitative predictions based on a 
toy model by Dekel & Burkert (2013). In the latter study, about 
half of the star-forming disk galaxies at z ~ 2 are predicted to 
evolve into compact star-forming “blue nuggets” due to violent 
disk instabilities before they are quenched into compact QGs 
(“red nuggets”). An observed population of “blue nuggets” has 
been proposed by Barro et al. (2013a, Barro et al. 2013b) to 
represent an evolutionary link, originating from extended disk 
galaxies and evolving into compact quiescent systems. 

This does not refute that galaxies also undergo further 
structural evolution after they are quenched. At least part of 
the size growth (Cassata et al. 2013; van der Wei et al. 2014) 
and evolution toward rounder axial ratios (Chang et al. 2013) 
has been attributed to (minor and/or major) dry mergers, 
and it is conceivable that similar processes contribute to the 
observed trend of increasing B/T toward the massive end for 
the quiescent population. 

6.2. AGN as the Driver of Quenching? 

We have demonstrated that the bulge mass of a system is well 
correlated with its quenched state and has a stronger predictive 
power of quiescence than the total stellar mass. The observed 
trends of / quen ch with total stellar mass, bulge mass, and disk 
mass as viewed among galaxies in different B/T bins are in 
good qualitative agreement with predictions from the Somerville 
et al. SAM. In the context of this model, the growth of the central 
supermassive black hole, which is the primary quenching agent 
for massive galaxies in this SAM, is tightly coupled with the 
growth of bulges through both merging and disk instabilities. 

If a black hole — bulge scaling relation is in place during the 
peak of cosmic star formation as it is in the present-day universe 
(Haring & Rix 2004), our observational results together with the 
model comparison could therefore hint at the bulge not being 
the causal link to quenching, but rather the most accessible 
observational proxy for the AGN acting as the quenching 
agent. 2 ' In detail, however, there are quantitative differences 
between the SAM and our observations, most severely in the 
highest redshift bin (1.5 < z < 2.5), where the observed 
quenched fraction exceeds the value predicted by the SAM by 
a factor of ~3.5. The latter difference could hint at a need 
for more frequent, efficient, or lasting quenching, a possible 
mechanism that we speculate about below. We also note that the 
same behavior is not necessarily a generic feature to all SAMs 
(see Appendix C). 

It is tempting to draw connections between the emerging 
bulges in massive MS galaxies out to z ~ 2.5 revealed by 
our analysis, and recent observational results based on deep 
adaptive-optics-assisted integral field data sets and grism spec- 
troscopy over the same redshift range. Forster Schreiber et al. 
(2014) found a high prevalence of powerful nuclear outflows 
in log (M) >11 galaxies driven by AGN, which appear to be 
absent in galaxies at lower masses. Along with star formation 
driven winds in the outer parts of the galaxies, such outflows 
could efficiently remove gas from the galaxies and, in this way, 
contribute to the quenching process. 


23 We note that Rosario et al. (2013) find X-ray signatures of AGN activity at 
these high redshifts to be most prominent among the star-forming population, 
most notably at the high-mass end, precisely where we see an upturn in the 
bulge fraction among SFGs. 
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Meanwhile, the 3D-HST and CANDELS legacy programs 
have yielded evidence for nuclear depressions in the Ha equiv- 
alent width in z ~ 1 SFGs (Nelson et al. 2012, 2013; Wuyts 
et al. 2013). At the highest stellar surface mass densities, star 
formation no longer appears to proceed in lockstep with the 
assembled stellar mass. Likewise, Genzel et al. (2014) report 
on ring-shaped Ha distributions in z ~ 2 SFGs, surrounding a 
more quiescent center where the dynamically inferred Toomre 
Q parameter significantly exceeds unity, owing to the emergence 
of a stellar bulge. Thus the Toomre stability criterion is satisfied 
in the central galaxy regions, which consequently could prevent 
the gas reservoir, if present there, from fragmenting and form- 
ing stars. While this result suggests that some causal connection 
between bulge growth and quenching may be at play, it should 
be noted (as is done also by Genzel et al. 2014) that gravitational 
quenching by itself does not expel the gas present, or stop the 
accumulation of a larger gas reservoir by continuing cosmologi- 
cal accretion. Additional maintenance modes might be required 
for a long-term shut down of further gas supply. 

7. CONCLUSIONS 

We analyzed the structural properties of a sample of 6764 
massive (>10 10 M 0 ) galaxies in the redshift range 0.5 < z < 
2.5, by exploiting the multi-wavelength CANDELS HST imag- 
ing data set in all five CANDELS/3D-HST fields. We carried 
out single-component (Sersic) fits and two-component (bulge + 
disk) decompositions, on stellar mass maps reconstructed from 
a resolved panchromatic SED modeling technique (Wuyts et al. 
2012, 2013), as well as on images of the //-band surface bright- 
ness distribution. In addition, we compared our findings to pre- 
dictions by the state-of-the-art SAM from Somerville et al. 
(2008, 2012; with extensions including disk instabilities pre- 
sented by L. A. Porter et al. 2014, in preparation). Our main 
results are the following. 

1 . At fixed stellar mass, QGs have overall higher Sersic indices 
and B/T ratios than SFGs as measured on their mass 
maps, in line with previous findings using monochromatic 
observations. We find an increase of Sersic indices among 
SFGs with increasing total stellar mass, with the median 
mass profiles increasing from ( n ~ 1.3) at 10 10 M 0 
to n > 2 above 10 11 M 0 . Two-component bulge-disk 
decompositions confirm that the same rising trend is present 
when considering the median B/T ratio of SFGs, which 
is rising up to ~40%-50% above 10 11 M 0 . The same 
characteristic behavior is seen at z ~ 1 and z ~ 2. 

2. Quantifying the same trends on the H- band light profiles, 
rather than the mass maps, the Sersic indices and B/T 
fractions are overall lower for SFGs, confirming previous 
non-parametric measurements for a subset of our sample 
(Wuyts et al. 2012). The emergence of bulges above 
10 1 1 M e in SFGs appears to also be slightly less prominent 
when viewed in light, which is consistent with the steepest 
color gradients (blue disks with red central bulges) being 
found among massive SFGs. 

3. The likelihood of a galaxy being quenched, as traced by 
the fraction of QGs, is better correlated with the bulge 
mass than the total stellar mass, and further shows no 
appreciable correlation with the amount of stellar mass in 
the disk component. The quenched fraction at redshift 1 is 
on average higher by a factor ~2 than at redshift 2. 

4. At a given total stellar mass, the quenched fraction exhibits 
a strong positive correlation with B/T , while different B/T 


bins are confined to a significantly tighter locus in a diagram 
of /quench versus Mb u i ge . These findings imply that the bulge 
mass of a system is the single observable parameter with 
the most predictive power regarding its quenched state, 
although a somewhat tighter constraint on the probability 
of quiescence can be obtained if the total stellar mass is 
also known. The same trend is seen over the full redshift 
range probed, with the distinction that quenched fractions 
are lower at higher look-back times. 

5. We find a good qualitative agreement between the SAM by 
Somerville et al. SAM and our observational findings. Since 
bulge and black hole growth are tightly coupled in the SAM, 
the strong dependence of / qu ench on bulge mass follows 
rather naturally in this model. Our observational results can 
in the context of this model therefore be interpreted as the 
bulge being the closest observable proxy to the underlying 
agent of quenching: the black hole. Quantitatively, the 
largest discrepancy between the model and observations 
is found in the highest redshift bin (1.5 < z < 2.5), where 
the observed quenched fraction is larger by a factor of ~3.5 
than predicted by the SAM. We note that the same behavior 
is not necessarily a generic feature to all SAMs. 
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APPENDIX A 

THE MEANING OF A SERSIC INDEX MEASUREMENT 

The combination of our one- and two-component fits on the 
mass maps enables us to empirically examine how the Sersic 
index relates the amplitude of the bulge (as parameterized by 
B/T) on the one hand, and the relative size-ratio to the bulge 
and disk component (R e ,B / //.,£>) on the other hand. Figure 6 
illustrates that, while B/T shows a clear correlation with n , there 
is no unique one-to-one translation between the two. Instead, 
the best-fit n to a composite bulge+disk system additionally 
depends on the size ratio of the two components. In other words, 
a galaxy’s Sersic index could be increased by placing more 
material in the bulge, but also by growing the disk at a fixed 
bulge size. The observed CANDELS galaxies occupy a surface 
in this three-parameter space (n, B/T, R c b / R e , i> ) that is in 
good agreement with what would be anticipated from Sersic 
fits to idealized, noise-free bulge (n = 4) plus disk (n = 1) 
profiles (squares and curves in the left- and right-hand panels of 
Figure 6, respectively). Evidently, for systems with B/T close 
to 0 or 1, the size ratio of the two components is ill-constrained, 
as one of them contains barely any mass. The galaxies in our 
full 0.5 < z < 2.5 sample span the full range of B/T values, 
and are located predominantly around bulge-to-disk size ratios 
ofR e , B /R e , D ~0.2. 

APPENDIX B 

COMPARISON WITH MEASUREMENTS ON H BAND 

Here, we investigate how the results of our structural analysis 
change when conducting the measurements on the H - band light 
images, rather than on the stellar mass maps. To this end, we 
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Figure 6. Left: size ratio between the bulge and disk component (R e ,B / Re,D) vs. B/T for the whole galaxy sample as measured in mass. The color indicates the 
measured Sersic index. Quiescent galaxies are shown as a plus (+), SFGs are shown as dots. Right: measured B/T values against Sersic indices as seen in mass. The 
histograms show the respective measured distributions of B/T , R e ,B/Re,D , and n among the sample for quiescent galaxies (solid black areas) and SFGs (hatched 
areas). The relationship between the Sersic index and B/T is not unique, but rather crucially depends on R e ,B/Re,D • The squares (left) and curved lines (right) 
illustrate the relation between {B/T, R e B / Re,D ) and n as inferred from noise-free idealized bulge + disk toy models, which are consistent with the trends seen for our 
observed galaxy samples. 

(A color version of this figure is available in the online journal.) 





10.0 10.5 11.0 9 10 11 9 10 11 

LOG (M.[M 0 ]) LOG (M Bulge [M Q ]) LOG (M D[sk [M 0 ]) 


Figure 7. Equivalent of Figures 1 (left) and 3 (right) using the results of structural measurements on the H- band light images. 
(A color version of this figure is available in the online journal.) 


present in Figure 7 the same figures as discussed in Section 4, 
now using the Sersic index and B/T values as inferred from the 
CANDELS H - band imaging. Likewise, we compute the bulge 
and disk mass as (B / T)hM* and (1 —(B/ T)h )M*, respectively. 

Overall, our analysis reveals a qualitatively similar mass 
dependence of n and B/T, and distinction between SFGs and 


QGs as inferred from the mass maps. In detail, however, modest 
changes in n and B/T are notable. While the median z ~ 1 
(2) SFG has (B/T)h < 0.25 (0.20), the typical bulge fractions 
increase to above 20%, and reach up to ~40%-50% at the 
massive end, once spatial M/L variations are corrected for. 
Likewise, the corresponding Sersic indices measured on the 
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Figure 8. Fraction of massive galaxies ( M > 10 10 Mq) that are quiescent in the SAM of Guo et al. (2013), as a function of total stellar mass, stellar mass in the bulge 
and disk separately, and as a function of halo and black hole mass. The bottom two rows show the results split by bulge-to-total ratio, for the z ~ 1 and z ~ 2 bins, 
respectively. 

(A color version of this figure is available in the online journal.) 


H - band imaging for SFGs below log(M) = 10.8 are consistent 
with the exponential disk profiles (see also Wuyts et al. 2011), 
but are slightly cuspier, as quantified on mass maps. This is in 
line with findings based on smaller subsets of CANDELS data 
by Wuyts et al. (2012) and Guo et al. (2012). 

The central 50th percentile intervals marked by the red and 
blue polygons are somewhat less confined in the plots based on 
stellar mass maps compared to the //-band results. We interpret 
this to be due to an additional source of random uncertainty 
introduced by the resolved stellar population modeling. The 
resolved stellar population modeling itself was motivated by 
the need to reduce the systematic biases associated with spatial 
M/L ratio variations. 

In Figure 7, the qualitative trends of / quenc h with total stellar 
mass, bulge, and disk mass are very similar. However, the lower 
two rows of Figure 7 compared to Figure 3 show that the bins 
of the lowest and highest B/T are more separated from each 
other in / que nch in light than in mass. This observation is also 
in line with the disks of SFGs being dominating by a younger 
stellar population than that of the bulge, shifting SFGs to lower 
B/T . The middle panels of Figure 7 show a larger scatter than 
the corresponding panels in Figure 3. 

APPENDIX C 

COMPARISON TO GUO ET AL. SAM 

In order to evaluate the extent to which the characteristic 
trends of the Somerville et al. SAM, as presented in Section 5, 
are generic to all SAMs, we considered an independent SAM by 
Guo et al. (2013). The Guo et al. (2013) model is rooted in the 
dark matter backbone of the Millennium Simulation (Springel 
et al. 2005), and its output tables with galaxy properties for 


snapshots of different look-back times are publicly available 
on the Virgo Millennium Database (Lemson & the Virgo 
Consortium 2006). We present the equivalent plots of Figures 4 
and 5 in Figure 8. As in the Somerville model, quenched 
fractions are rising with increasing mass, increasing bulge mass, 
and decreasing redshift, but when broken down in bins of B/T, 
significant differences are notable. Specifically, / que nch does not 
monotonically increase with increasing B/T at a given stellar 
mass. 
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